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Abstract 

We attempt to build systematically the low-energy effective Lagrangian for the Einstein-Cartan 
formulation of gravity theory that generally includes the torsion field. We list all invariant action 



Q_i terms in certain given order; some of the invariants are new. We show that in the leading order 

CD 

the fermion action with torsion possesses additional {7(1)l x f (1)r gauge symmetry, with 4+4 

t^J- ' components of the torsion (out of the general 24) playing the role of Abelian gauge bosons. The 

>' 

bosonic action quadratic in torsion gives masses to those gauge bosons. Integrating out torsion one 
. obtains a point-like 4- fermion action of a general form containing vector- vector, axial- vector and 

t^J- \ axial-axial interactions. We present a quantum field-theoretic method to average the 4-fermion 



interaction over the fermion medium, and perform the explicit averaging for free fermions with 
given chemical potential and temperature. The result is different from that following from the 



^ \ "spin fluid" approach used previously. On the whole, we arrive to rather pessimistic conclusions 

H ■ 

on the possibility to observe effects of the torsion-induced 4-fermion interaction, although under 
certain circumstances it may have cosmological consequences. 
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I. INTRODUCTION 



Recently, there has been some renewed interest in torsion appearing in the Einstein- 
Cartan formulation of General Relativity jl-6|, and in physical effects it may imply, see 
e.g. Refs. A general drawback of these interesting studies was certain lack of 

systematics, in particular not all a priori possible invariants of a given order containing 
torsion were considered. The first aim of this paper is to fill in this gap and to introduce 
torsion including its interaction with fermions in a systematic way, in the spirit of the 
low-energy derivative expansion Q. In the .owest order of this expansion, torsion induees 
a local 4-fermion interaction [1| which may affect cosmological evolution at very early 
times and high matter density |lll-ll4j. or be detectable from high-precision data at later 
times. Therefore, the second aim of this paper is to evaluate the contribution of the general 
torsion-induced 4-fermion interaction derived in the first part, to the stress-energy tensor 
for possible applications. We think that such evaluation in the past, based on "spin fluid" 
ideas, has been unsound. 



General Relativity with fermions is a theory invariant under i) general coordinate trans- 
formations (diffeomorphisms) and it) point- dependent (local) Lorentz transformations. The 
standard way one introduces fermions is via the Fock-Weyl action 
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S { = i Jd 4 xdet(e)^(Ve A ^ A V^-V^e A » lA y), V, = d, + \uf[ lBlc ], (1) 

where $ is the 4-component fermion field assumed to be a world scalar, 7^ are the four 
Dirac matrices, u^ c is the gauge field of the local Lorentz group, called spin connection, 
and e A ^ is the contravariant (inverse) frame field. In fact there are other fermion actions 
invariant under i and ii, to be discussed below, but Eq. (CQ) is generic. 

To incorporate fermions, one needs, therefore, the gauge field and the frame field e M , 
which are a priori independent. Therefore, the bosonic part of the General Relativity action 
must be also constructed from these fields. We are thus bound to the Einstein-Cartan 
formulation of General Relativity. 

In this formulation, the lowest-derivative terms in the bosonic part of the action invariant 
under i and ii are 

Mp r 4 / , 1 A -[-AB CD 1 nXav t-AB A B , \ / Q \ 

WirJ V ~ 4 ^abcd J^k\ e^v - 2 e e^e u + ...j (2) 



where D AB = (d^ + u^) AB is the covariant derivative, and T AB = [D^D U ] AB is the curvature. 
The first term is the cosmological term, det(e) = (l/4!)e KA/ " / €abcd e^ef e^e B = yf—g, 
with A = 5.5 • 10~ 84 Gev 2 , the second term is the would-be Einstein-Hilbert action with 
Mp = 1/ VG = 1.22-10 19 GeV being the Planck mass, and the third term is the P- and T-odd 
action first suggested in Refs. 17], Il8] with iota being so far a free dimensionless parameter. 
In the context of canonical gravity the inverse of t is sometimes called the Barbero-Immirzi 

called the "Hoist action" [l9] although in fact they 



171 . Il8j . In a quantum theory one writes the amplitude as 



parameter, and the third term in ([2] 
were introduced much earlier 
exp(iS). 

Eq. 02]) is quadratic in u^, therefore the saddle-point integration in uj^ in the assumed 
path integral over u^, e„, is exact. We define the antisymmetric torsion tensor as 



rpA A 



\ ((D,e u ) A - {D u e,) A ) . (3) 



Under Lorentz transformations it behaves as a 4-vector, and under diffeomorphisms it be- 
haves as a rank-2 tensor. In the absence of fermions, the saddle-point equation arising 
from the first variation of Eq. (jSJ) in 24 independent variables u AB results in 24 dynamical 
equations T A U = 0. This set of equations is linear in u^, and its solution is the well known 

rf B (e) = \e AK {d,e B - d K e B ) - \e B %d»e A - d K e A ) - \e A ^e%d K e c x - d x e c K ). (4) 

Quantities with a bar refer, here and below, to this zero-torsion case. [If two flat in- 
dices A,B,C, ... appear both as subscripts or both as superscripts we sum over them with 
Minkowski signature t] = (1, — 1, — 1, — 1); e 0123 = 1, £0123 = — 1.] It is well known that in 
the zero-torsion case the second term in Eq. fl5]) reduces to the standard Einstein-Hilbert 
action —yf—gR where R is the standard scalar curvature built from the Christoffel symbol 
r re v = \{dn9\n + d >~9Kfi - d^g K \), and the third term reduces to -|-R K A,^e KA ^ = 0. 

With fermions switched in, the torsion is nonzero even at the saddle point since from 
varying Eq. ([[]) one obtains T ~ J, where J ~ is a fermion bilinear current. Neither is 
it generally speaking zero, if terms with higher derivatives are added to the bosonic action, 
see Section VI. C. 

Apart from full derivatives, the deviation of the second and third term in Eq. ($Z§ from its 
zero-torsion limit is 0(T 2 ), and the deviation of the fermionic action ([1]) from its zero-torsion 
limit is 0(T), see Eqs. fl60|61|l55|) below. However, these terms are not the only ones that 



can be constructed from the requirements i and ii, and there are no a priori reasons why 
other terms should be ignored. 

The minimal actions (TTP|) are at best the low-energy limit of an effective theory whose 
microscopic origin is still under debate, since the action is non-renormalizable and, worse, 
non-positive definite. Therefore, the best we can do in the absence of a well-defined quantum 
theory is to write down a derivative expansion for an effective low-energy action satisfying 
the requirements i and ii with arbitrary constants, to be in principle determined or at least 
restricted from observation or experiment. A future fundamental, microscopic theory of 
gravity will be able to fix those constants. Unless proved otherwise, the derivative expansion 
in the effective action is assumed to be a Taylor series in d 2 /Mp in the bosonic sector; in 
the fermionic sector odd powers of d/Mp are also allowed. 

For completeness, we include in our consideration invariants that are odd under P, T 
inversion. Since these discrete symmetries are not preserved by weak interactions, and the 
effective low-energy gravity may imply integrating out high-momenta fermions, we do not 
see the principle why such invariants should be avoided. 

II. ORDER OF MAGNITUDE ANALYSIS 

In the Einstein-Cartan formulation of gravity, which as we stress is unavoidable if we 
wish to incorporate fermions, the frame field and the spin connection are a priori 
independent variables. One expands the action in the derivatives of e M and u^, preserving 
the diffeomorphism invariance and gauge invariance under local Lorentz transformations, 
the only two symmetries requested. The spin connection is a gauge field that transforms 
inhomogeneously under local Lorentz transformations, hence it can appear only inside co- 
variant derivatives. 

At the saddle point (in other wording from the equation of motion) the spin connection 
is expressed through the frame field, = ui^ ~ e~ x de, see Eq. (@|. Suppressing the 
indices and omitting the frame fields assumed to be of the order of unity, one can present 
symbolically the general as 

lo = lo + T, T~<9e, (5) 
where T is the torsion field fl3]); the precise relation is given by Eq. fl35|) below. From the 
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point of view of derivative counting, T, uj and hence u itself are all one-derivative quantities. 
From the point of view of the gauge group, u transforms inhomogeneously and hence must 
always come inside covariant derivatives, whereas T transforms homogeneously and hence 
terms of the type T 2 and the like are allowed by gauge symmetry. 

Let us classify the possible action terms. There is only one zero-derivative term, the 
invariant volume or the cosmological term. One-derivative terms are absent in the bosonic 
sector. They appear only in the fermionic sector from the Dirac-Fock-Weyl action which 
we fully analyze in Section IV: there is a unique term of the type ^>(d + (D)^ and four terms 
of the type ^T*. 

There are precisely seven two-derivative terms: two terms linear in the curvature J 7 , 
presented in Eq. 02]), and five terms quadratic in torsion, fully listed in Section VI. A. Since 
by definition T ~ Ve, all terms quadratic in torsion are quadratic in the derivatives. There 
are no other terms quadratic in the derivatives in the bosonic sector, and this is an exact 
statement of this paper. 

Turning to four-derivative terms, there are in general terms of the type J 7 J 7 ~ R 2 , T(VJ-), 
(VT) 2 , T 2 (VT), T 2 T and T 4 . 

Omitting the cosmological term we write down symbolically the effective Lagrangian as 

£ eff = M$R + Vd% + *T^ + M 2 P T 2 (6) 
+ R 2 + T(Vi?) + (VT) 2 + T 2 (VT) + T 2 R + T 4 + 0(1/ M 2 ) 

where the Lagrangian for matter is represented by the fermionic source. The equations of 
motion are obtained by varying the action ([6]) with respect to the fields involved. Assuming 
h = c = 1 and the dimensionless metric tensor g^ v , the dimensions of the quantities in 
Eq. 06]) are the usual R ~ 1/cm 2 , T ~ 1/cm, ~ 1/cm 3 . In the leading order one gets 
from the first two terms the standard estimate for the curvature generated by matter, 

- < ggg > P 4 ( , 

Mp ~ M 2 ' 1 ' 

where p is the characteristic momentum of matter, be it temperature, cubic root of density, 
or mass - in fact the largest of these. Being substituted back into the action, this estimate 
shows that the first two terms of the action (J6j) are of the order of p 4 , and that the R 2 term 
is a tiny p 8 /M p correction. 
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Excluding torsion T from the 3 d and 4 th term in (J6]) gives 

< ^ > p 3 



(8) 



< > p 



M 2 M 2 p 

This is the leading post-Einstein correction, and we analyze its most general structure and 
its effect in this paper. From the estimate (jSJ) we see that other terms in Eq. ([6]) give even 
smaller corrections: 

T 2 ~ T{VF) ~ (VT) 2 ~ (9) 



THWT) ~ T J ~ — T ~ — — 

Nevertheless, we list for completeness all 10 possible terms of the type T 2 in Section VLB 
and all 4 possible terms of the type T(VJ r ) in Section VI. C. 

III. GENERAL FRAMEWORK 

In this section we introduce the basic variables and make sign and normalization con- 
ventions. To simplify the algebra, we temporarily deal with the Euclidean signature where 
the Lorentz group SO (4) acting on flat indices A,B,C,... is locally isomorphic to the direct 
product S77(2)l X S77(2)r.. We return to Minkowski signature in the final results. 

A. SU (2)l x SU(2)r subgroups of the Lorentz group 



The 4-component Dirac bi-spinor field in the spinor basis is 

/ 

q? = 



* f = te,xi). (io) 



4 

The Lorentz SU{2)\ J x S77(2)r transformation rotates the 2-component Weyl spinors: 

r -> (u L );^, £ -> 4 (ui)[, u L g su(2) L , 

x a -> (^)iv^, xl xj (4)^ % e 5C/(2) R . (li) 

In the spinor basis the (Euclidean) Dirac matrices are 

/ 

ai = {l 2 ,±ir i ), [a\y = a A . (12) 



a~ A 

lA = | | , 75 

o-i o 



A u / V 



i 2 o 

-1, 



where r % (i = 1,2, 3) are the three Pauli matrices. We introduce the commutators 

v-H — d_ 1 ( + — + -\ _ i i ( y-H — \t v-H — 

^AB — 2 \ a A (J B ~ <7 B°'A J — ~VAB T j \ AB ) ~ ABi 

^ab = \ (oa<J% - (Jb^a) = -V 1 abt\ (s A ^) f = (13) 

where 77,77 are 't Hooft symbols. They are projectors of so(4) to the two su(2) subalgebras. 
The basic relations for 't Hooft symbols are 

VabVcd — 5ac$bd — Sad^bc + cabcd, 
VabVcd = $ac8bd — 8ad8bc — £abcd, 

v\bVab=^ 1 \ v\bVab = ^ 1 \ rf AB ff AB = Q. (14) 

With Euclidean signature, there is no distinction between upper and lower flat (capital 
Latin) indices, in particular, e 1234 = €1234 = 1. 

The following commutation relations are helpful and will be used below: 

^BC a A ~ Va^bC = ~ 2i (5BAVC ~ ^CA^b) , 

^BC a A ~ a A^BC = ~ 2i ($BA(Tc ~ ^CA^b) ■ ( 15 ) 

B. Tetrad 

We introduce the frame field in the matrix form: 

4 = 4°% (i6) 

Under Lorentz transformations, the frame field transforms as 

e+ URe^Ul e~ U L e-Ul (17) 

Under the general differentiable change of the coordinate system x^ — > x >fl (x) the frame field 
transforms as a world vector: 

e±(x)^ef(x'(x))—. (18) 
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C. Covariant derivatives and curvatures 



The requirement that the theory is invariant under local (point-dependent) SO (4) Lorentz 
transformation demands that a 'compensating' gauge field ui AB = —u BA , called spin con- 
nection, must be introduced. With its help, one constructs the covariant derivatives, 

Df = d,8 AB + u AB , (19) 
V; + = d,l 2 --u AB H- A i = d^U-iL*^ = (20) 



v 2 

i 

■ fi —/in -ft —£ n 2 



V+- = U - -A B ^\b = d^-iRl- = d,-iR^ (21) 



where L(R) are left (right) connections, 



^ B = -\ L ^ab-\rWab- (22) 



Inversely, 



In matrix notations 



4 = -\v\ B "f, R; = -\rf AB uf. (23) 

W = L^-, R, = Kj-. (24) 
Under gauge transformation ( fTTj) the covariant derivatives transform as follows: 

V-+ ->• U L V~ + Ul V + ~ ->• U R V+-Ul (25) 

corresponding to the usual gauge transformation of the connections, 

L„ -> U L {L^ + i5 M )c4, -> U R {R P + id^Ul (26) 

The commutators of the covariant derivatives are curvatures: 

(27) 
(28) 

(29) 

The SO (4) curvature is decomposed accordingly into two pieces transforming as the (3, 1) © 
(1,3) representation of the SU{2)\ J x S77(2)r group: 

= ~\ F U L ) Vab - \ FUR) Vab (30) 



[dm** 


-pAB 
•'fin 


- f) f , AB 


r) i , AB 


, , .AC. ,CB 






= F~ + 

flU 


= {d^Lv - 


d y L^ - 


i[L^L u }) = 




*[v+-vri 


= F + - 

flU 


= {d^R u - 


duRfj, — 


i[R P Ru\) = 


-f; u (R)j- 
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where 



f; u (l) = o^-d^ + ^u^ 

F; U {R) = d,Rl - + <«*J£J2* (31) 

are the usual Yang-Mills field strengths of the SU(2) Yang-Mills potentials and R l . 
These field strengths are projections of the full curvature: 

F+~(R) = F; U {R) T - = (32) 

D. Torsion 

The antisymmetric combinations 

(V M + e; - e;V+-) - {V- + e- - e- V+") = (pfg - D* B ef) a A = 1T*<T A , 
(V+"4 - - (V+"eJ - <V;+) = - Dt B e*) a\ = 2T>+ (33) 

define the torsion field T A U = | ((D^e^)" 4 — (Z^e^m = (D^e^)" 4 . It is a 4- vector with 
respect to Lorentz transformations, and an antisymmetric rank-2 tensor with respect to 
diffeomorphisms. Contracting it with e Aa one gets a rank-3 tensor T^ x = —T Ufl x = T A e . 

The torsion tensor T^ v x has 24 independent components. It is convenient to decompose 
T^ u x into the totally antisymmetric part related to an axial vector (a K , 4 components), the 
trace part related to a vector (v K , 4 components), and the rest 16 components (t A ) subject 



to constraints 



2l|: 



= \ det(e) a K e Kflvp g pX + ^v K 5^S x u] + j^, (34) 



where, inversely, 

rp rp f_l i A T'-^-L T 1 ^ n^P T 

' 1 a/3,75 U K — 1 Kfl 5 >K — 1 ~ 1 ~°[fM 1 u]p 9 1 p\w] 



4det(e) 

The reduced torsion tensor t ,f satisfies 4 constraints e KAMV t^t- = and 4 constraints t Mp p 
0, therefore it has 16 degrees of freedom, as it should. 
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The general 24-component spin connection can be presented as a sum of the zero-torsion 
part (j3J) and the torsion part: 



, ,ab _ -AB 



+ (T^ + T%, + V) (35) 



,- AB , 2 f ABCD p C p I>„K , ^ JB„A]k„, , 2 „Aa /j. , f , f \ 



E. Affine connection 

In the Einstein-Cartan formulation, the tetrad e A and the spin connection ui AB are pri- 
mary fields, whereas the metric tensor g^ u and the general affine connection V* are sec- 
ondary, defined through the first pair. 

The quantity D AB e B = (D fl e u ) A is a vector in the flat space, therefore it can be decom- 
posed in the frame field e A which forms a basis in the flat space. We denote the expansion 
coefficients by T* u , 

(D,e v ) A 1 J%e}, (36) 

which serves as a definition of the affine connection It is equal to the sum of the standard 
Christoffel symbol (or Levi-Civita connection) to which the general affine connection reduces 
in the zero-torsion limit, 

f £„ = \g XK (d,g KV + d u g K , - d K g, u ) , (37) 

and the torsion part, 

r;, = r^ + T^ + T^ + T^ A . (38) 

F. Affine curvature and Riemann tensor 

Eq. (|36|) can be rewritten as 

(V,T x e AX = -to AB e B «, (39) 

where 

(v^ = d^ x + r; x (40) 
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is the standard affine covariant derivative. The commutator of two covariant derivatives 
defines the affine curvature tensor 

NMl = R\„ = o,v: x - d u r; x + r; p rt x - K p r; x , (4i) 

which is related to the curvature (127|) built from the spin connection: 

R\,u = e AK e B x T£. (42) 

The all-indices-down curvature tensor and the generalization of the Ricci tensor are 

RkX^v — e£ ef T AB , R K \ — R Kf i,Xu 9^ = e Bv T x ® . (43) 

Let us denote with the bar the curvature defined by Eq. (14 ip but built from the symmetric 
Christoffel symbol f (l37j) . R K \^ U is then the standard, zero-torsion Riemann tensor satisfying 
the following relations: 

RkX,/ii/ R[mv,kX i 

RkX,/iu Rk[i,vX Rku,Xij, , 

RkX,^ £ KXfJ/U = 0, R K \ = R\ K . (44) 

These relations are not valid in the general case for the un-barred curvature R K x,^ ( 143]) if 
torsion is nonzero. 



IV. THE FERMIONIC ACTION 

In this section, we construct all possible bilinear fermion actions with zero and one co- 
variant derivatives. To make sure that we do not miss any terms, we prefer to use the 
two-component formalism, see Section III. A. Any action is, in principle, allowed that is i) 
diffeomorphism- invariant and ii) invariant under local Lorentz transformations (fTTj) . The 
first requirement means that, if only covariant (lower) indices are used for e^, V^, , 
they must be all contracted with the antisymmetric e in order to compensate for the 
change of coordinates in the volume element d A x. The second requirement means that, be- 
cause of the gauge transformation laws ( TTT|I17I)25[) , one has to alternate subscripts 'plus' and 
'minus' in the chain. Fermion operators are at the ends of the chain, such that the 'plus' 
is followed by the Weyl field ip whereas the 'minus' is followed by the Weyl field x- On the 
contrary, ip* is always followed by a 'minus' whereas x^ is followed by a 'plus'. 
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A. Zero-derivative terms 



The zero-derivative fermion action can have only two structures: 



si = T[ fd 



xie 



d xi det(e) 



(45) 



Both terms are Hermitian if / i an d /02 are rea l ; since ip, ^ and X-> X^ are Grassmann 
variables. The two terms transform into one another under parity transformation. Therefore, 
if parity is not broken /01 — — /02 — m - In this case the two terms combine into the mass 
term: 

Sl=i J d 4 x det(e) m * f *. (46) 



B. One-derivative terms 



Since all covariant indices must be contracted with e KXflu , the total number of covariant 
indices belonging to e^,V^ T must be four. In one- derivative terms, one covariant index 
belongs to V^ T , therefore the number of frame fields must be three. It means, first of all, 
that the number of 'pluses' and 'minuses' is odd, therefore all possible one- derivative terms 
are necessarily off-diagonal in the Weyl fields ip,X- It is sufficient to consider operators of the 
type (x ■ ■ ■ VO since the opposite order (ip^ ■ ■ - x) wm b e obtained by hermitian conjugation. 

A priori one can construct many terms satisfying these constraints, however all of them 
can be reduced, using the algebra from Section III. A, to the following four terms: 

eK A^ (x t e + V -+ e - e + V , )j e ^( % t V +- e + e - e + V , ) _ (4?) 

The covariant derivative in these strings acts either on the Weyl field or on the tetrad. 
One can commute the derivative to the utmost right position, therefore the only term with 
the derivative of the fermion field is the usual Dirac term, e KXllu (x^ e~^'V~ + ip). All other 
terms contain derivatives of the tetrad, which should be antisymmetrized in the covariant 
indices to make world tensors under diffeomorphisms. This forms the torsion tensor (}33l) . 
Therefore, apart from the Dirac term, Eq. (j4T|) describes the following three structures: 

e KX ^ tfT+e-etl*), e^^e+T^), e^x^T^) (48) 
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where = T^a^- In addition to the four terms in Eq. (H7|) one can consider 
^(x f e»Tr(T+ e ;) and e^^-ty T*(I£e+). 



These two terms are in fact identical and equal to I e KA /^ n^e+e^T^^) — (x^T^e^ e\ ■ 
belonging to the set fj48|) . Further on, the three terms in Eq. (j4*8]) are not independent as 
there is an algebraic identity 



leaving us with only two terms with the torsion field, say, the first and the last term in (|48|) . 

In principle, one can build diffeomorphism-invariant actions using two Levi-Civita sym- 
bols e KXfMU contracted with one V^ T and seven but divided by det(e) such that the expres- 
sion is again invariant under the change of coordinates, together with the volume element. 
However, all such expressions are in fact identical to linear combinations of invariants listed 
in gTD. 

The resulting three independent fermion actions can be presented in a more simple form. 
We use the decomposition of the torsion tensor ( 1341) and notice that actually only 8 out of 
the possible 24 components of the torsion field couple to fermions in this order: the traceless 
symmetric part £ * decouples. Indeed, one has: 



eK A^ (x t e + e - e +v-+V>) = -6det(e)( X t e +^V; + V0, 
e^ u (^T+e-ei^) = -8det(e) + a M ) (x f e + ^), 

e ^ (x t e + e - T + V , ) = _8det(e) - a^j (x^). 

Therefore, the most general one-derivative fermion action has the form 

S[ = [ d 4 x det(e) \fio(x f e +X V^) + f u a A (x + e+V) + fu ^(x t e +A V)l + h.c, (49) 



where fu are arbitrary complex numbers. 

It should be noted that the operator in the first term contains the full spin connection 
which, according to Eq. ( I35l) . can be written as a sum of the zero-torsion part ui^ (@J and 

terms proportional to torsion. The difference of two operators is from ( |35|34|) identically 

e ±A (V^ - Vf=) = e ±x (-v x T a A ) (50) 
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where Vj 1 * 1 is the covariant derivative computed in the zero-torsion limit. We have 

S[ = J d 4 x det(e) [^o(x t e +A V^ + ^) + g x a A (x f e + V) + 92 ^(x f e + V)] + h.c. , (51) 
where the new couplings are #0 = f w , g x = f u - fio, 92 = /12 - fw- 

C. Hermitian action 

We now add explicitly the Hermitian conjugate action in Eq. (|51|). It is straightforward 
for the second and third terms; changing the order of the fermion operators brings the minus 
sign. The Hermitian conjugation of the first term is more involved. We first take its complex 
conjugate, interchange the order of fermion operators (%, 1/)*), integrate by parts and get 



det(e) (Ve +A V A + ^) h c ^ = (^V^+e^ det(e) X 



We next write e~ A det(e) = — (l/6)e Aa/?7 e~e~^e~ and drag to the right through this 
expression. Since torsion in is by construction zero, it commutes with the tetrads 

owing to Eq. (!33j) . and we get 

tV^ + e- A det(e)x) = (Ve- A det(e) Vj"x) • 



We, thus, obtain an explicitly Hermitian fermion action written in the 2-component spinor 
Weyl form 

S[ = J d'x det(e) [go(x j e + ^- + ^) + ^(V^V^x) (52) 
+ a, ( 9l ( X j e + ^) - gl^e-^x)) + (g 2 (x t e + ^) - g^e^x)) 



The constant go can be made real by redefining the overall phases of if), if)' and XiX ■ 
Indeed, if the argument of go is a (go = \go\e ia ), the phase rotation if) — > exp(—ia/2)ip, ifr — > 
ify exp(ia/2), x exp(ia/2)x, x 1 ~ ^ X*exp(— iot/2) obviously makes go real, and it can be 
further on put to unity by rescaling of the if), x fields. Therefore, we can put go = 1 to make 
the Dirac kinetic energy term standard. 

Finally, recalling the definition of the bi-spinors f llOp and the Dirac matrices (1121) we 
rewrite the action f[5"3"j) in the 4-component Dirac form: 

S{ = J d A x det(e) *Y + a„(g x + 2i + 7s) + + 0+75)] * (53) 
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where 



± _ 01,2 ± 9*1,2 
#1,2 _ o : 



r 



e Afl 7A, 



V f , = d fl l + l[~f A -f B ]u 



AB 



(54) 



The action f[53|) is by construction and manifestly Hermitian 



In the 'minimal model' often discussed in the literature jl, 5, 22], the only source of the 
fermion interaction with torsion is the Dirac term ([T]) with the full spin connection including 
its torsion part. In this case gf = —1, and all other constants are zero, therefore only 
the axial part of the torsion couples to fermions. The term proportional to g% was first 
considered in Ref. |tJ although in another form, see also its discussion in Ref. [6|. These 
authors take the following Lagrangian generalizing the 'minimal model': 

1 



= Jd A x det(e) 



W^V^ — (PvP) V* 



2 ' ' p 2 

Integrating the second term by parts we bring this Lagrangian to our form (153]) with the 
particular values of the constants: gf — — 1, g% = a, the rest being zero. 

The complete list of four one-derivative fermion-torsion actions ( 1531) was presented in 
Ref. 23] where in addition nine terms with one extra derivative were suggested. 



D. Prom Euclidean to Minkowski signature 



The standard dictionary translating Euclidean into Minkowski variables (see, e.g. 24j ) 
reads: 

7° = 74B, i = -Hei y = ^E, * = i*t» 7 5 = S = lS E . 

Therefore, we obtain from Eqs.( )46|53|) the full Minkowski fermion action with zero and one 
derivative: 



S f = f d A x det(e) * [i^ (p„ + a p {g^ - ig+j 5 ) + - ^ 2 V)) - m 



(55) 



Assuming a M is an axial and t> M is a vector field, the terms with g 1 and break P-parity; 
other terms are parity-even. We remind that for a quantum amplitude one takes exp(iS). 

V. TORSION AS AN ABELIAN GAUGE FIELD 

If parity is conserved, one has to put g{ = g£ = in Eq. (155]) . We denote the two 
nonzero constants that are left as gf = g a (it is real by construction fl54|) ) and g% = —ig v 
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(it is purely imaginary), and rewrite the fermion Lagrangian with torsion as 

/- 1 + 7 5 1~7 5 \ 
#1" [Dp-Hp — ir M — - — I*, 1^ = g v Vp + g a a^ r M = (/„ u M - V ( 56 ) 

This Lagrangian is clearly invariant under the Abelian x £/(1)r gauge transformation 

# L = V^e iQ ^V, VR = X^e^ X , l^l^ + d^a, r^r^ + d^. (57) 



This invariance, in the U(l)y x U(1)a form, has been previously noticed in Ref. 25]. 

Therefore, the inclusion of torsion is equivalent to promoting the fermion part of the stan- 
dard gravity invariant under the Lorentz S77(2)l x S77(2)r gauge group to being invariant 
under the larger U(2) L x U(2) R gauge group! 

If this gauge symmetry is preserved by the bosonic part of the action, it has to depend 
only on the curvatures ff£y = dpl v — <9„Z M , fffi = d^r^ — d^r^. Linear terms are zero, so 
the expansion starts with quadratic terms in The vector part, v^, is then identical to 
the photon. It should be mentioned that, since it is an Abelian field, one is free to ascribe 
arbitrary coupling constants or 'charges' with which this field interacts with various fermion 
species. To be separated from the photon, this field has to have other 'charges' with respect 
to fermion species, and be massive. The only thing we know is that its mass must be larger 
than the experimental restriction on neutral intermediate bosons, of the order of 1 TeV. The 
same applies to the U(l) axial boson a M or their linear combinations Z M , r M . 

The appearance of mass terms for 1^, r M fields means breaking of the {7(1)l x /7(1)r gauge 
symmetry, which can be either explicit or spontaneous by some kind of a Higgs effect. As 
we shall see in the next section, adding terms quadratic in torsion implies explicit breaking 
of the {7(1)l x C/(1)r symmetry. In principle, there is nothing wrong about it as the gauge 
symmetry is Abelian. Probably, spontaneous breaking would be more aesthetic but in the 
absence of the microscopic theory we can only speculate about it. 

From the viewpoint that torsion fields a^v^ are just another set of gauge vector bosons 
interacting with fermions, we do not see compelling reasons why their masses should be 
of the order of the Planck mass, as suggested by the 'gravitational' approach to torsion: 
with our present lack of deeper understanding the masses can be anything beyond the 
phenomeno logically established limits (26 ]. 

If parity is not conserved (meaning g±, g£ are nonzero), one can still consider 1^ and as 
compensating gauge fields. However, then they have to be complex, and compensate point- 
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dependent real dilatations of the chiral fields ip, x an d not only their phases. A discussion 
of this interesting topic lies beyond the scope of the paper. 



VI. THE BOSONIC ACTION 

Taking a purely phenomenological stand, one may inquire what terms in the bosonic 
action can be written that preserve i) diffeomorphism-invariance and ii) invariance under the 
gauge Lorentz group. In this section, we give the full list of invariants quadratic in torsion, 
invariants quadratic in curvature, and invariants that are linear in torsion but containing 



A. Invariants quadratic in torsion 

A general way to construct quadratic invariants is to consider the following invariant 
under the diffeomorphism, 



and to contract the flat indices into a Lorentz-group scalar. Since torsion is antisymmetric 
in world indices, this expression is antisymmetric in [CD] and [EF], meaning that the pairs 
of the frame fields belong to the 6-dimensional representation of the SO (A) Lorentz group 
or to the (3, 1) © (1, 3) representation of the SU(2) L x SU(2) R group, while T A belongs to 
the (2, 2) representation of that group. 

The direct product of T A and e^e D ^ belongs to the (2, 2) <g> ((3, 1) © (1, 3)) = (2, 2)i © 
(2, 2) 2 © (2, 4) © (4, 2) representations, which should be multiplied by the same combination. 
There are 5 singlets arising from 



(2, 2)! © (2, 2) ls (2, 2)i © (2, 2) 2 , (2, 2) 2 © (2, 2) 2 , (2, 4) © (2, 4), (4, 2) © (4, 2). 



Therefore, there are precisely five linear independent invariants which we write as 



VR. 



K A[CD],B[EF\ = det(e) ( T A e ^ e D, )(T B £ Ea £ F^ 



(58) 



K X = det(e) T A T A gT tf* = det(e) (- 
K 2 = det(e) T A e A » e Ba g vf) = det(e) v> x v i 
A 3 = det(e) Tt T*« e Aa e B » = det(e) ( | . 
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^5 = fT* T£ T% e A K ef = det(e) | a% - ± t Q/3 A V, A e a ^. (59) 

In the last column we used the decomposition of the torsion tensor fl34l) . The last two 
terms are P, T-odd; the first three are even. The first four terms have been known for a long 
time as they emerge from the leading Einstein-Cartan terms (JSJ, see below. To the best of 
our knowledge, the fifth invariant appears for the first time in the very recent Ref. [27 ] . 

We now recall that in the Einstein-Cartan formulation there are two leading terms linear 
in the curvature, see the second and the third terms in Eq. ([2]). Following the general 
strategy we split them into a piece that survives in the zero-torsion limit, plus corrections 
from torsion. One has 0, 2s[ 29]: 

1 ^ e ABCD T AB e C e D = ^ R + Ki _ 4K2 + 2K3 + 4d ^^ v ,^ (6Q) 
1 ^ -pAB e A e B = ^ _ ^ ^ all j (61) 

We see thus that four out of possible five terms quadratic in torsion are induced by terms 
linear in curvature, with concrete coefficients. We shall, however, consider the general case 
where the terms i^i...s are included in the bosonic part of the action with arbitrary real 
coefficients fci...s: 

S2= fd 4 X^J2 k mK m (T) (62) 

= J d\ det (e) - L M 2 aa a% + 2M 2 av a% + M 2 VV v»v, + + - M e 2 tt t a/3 A V,a ^y) • 



The new constants having the meaning of the masses squared of the torsion fields are linear 
combinations of the original constants: 

M 2 _ ~ fci) M 2 _ 4(fc 4 + h) 2 _ 2fci + 3fc 2 + fc 3 
M2 = 2(2fcj + fc 3 ) M2 _ 2(2A; 4 - fc 6 ) 



9 ' '" 9 



The second and fifth terms are P,T-odd, the rest are even. The first three terms in Eq. ( l62l) 
are mass terms for the a^v^ bosons or for their linear combinations l^r^, that break 
explicitly the U{\) x £7(1) gauge symmetry discussed in Section V. 
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The system must be stable with respect to small low-momenta fluctuations of torsion 
about the flat space. It means that all eigenvalues of the mass matrix fl62|) must be positive. 
This condition requires that 

M 2 vv + M 2 aa >0, M 2 vv M 2 aa > Ml, Ml > Ml t (63) 

which is satisfied in a broad range of the constants h\^. See also the discussion of the 
positivity of the mass matrix in Ref. j^. 

In the 'minimal model' corresponding to extracting torsion terms from the leading-order 
action (J2J) only, see Eqs. fl60|61l) . one obtains 



2 hi = J_Ml M 2min = 8M| M 2m in= _8M| 

316tt' vv 316tt' u 916tt' 

M 2 min = § Mp M 2min = _ 8 Mp 

3 16vr' Ktt '9 16vr' l ° ] 

where the iota parameter i is the coefficient in front of the P, T-odd action (1611) . We coincide 
in this table of masses with Ref. 



301 ]. after adjusting the normalization. 



The eigenvalues of mass-squared matrix for a^, are ±vT+?(8/3)M|,/167r. A check of 
the above algebra is that at purely imaginary values i = ±i the eigenvalues are zero. Indeed, 
at these values the self-dual or anti-self-dual combination T AB ± %\e ABGD Tqd drops out of 
the action (T5]). 

At real values of the iota parameter one of the eigenvalues is always negative. It means 
that the path integral over a^, fields strictly speaking does not exist, therefore the 'minimal 
model' cannot be complete. 

B. Invariants quadratic in curvature 

Such terms arise from the diffeomorphism-invariant structure 

G ^B][GD]{EF][GH] = det(e) -pAB jrCD & E a ^ g G 7 & m ( q 5) 

belonging to the 6 ® 6 ® 6 <8> 6 representation of the Lorentz group, out of which one can 
extract 10 independent Lorentz-group invariants. Here is their list, expressed through the 
full Riemann tensor ([43]): 



Gl ~ 16deIM V 6 6 ^ e * H) ~ V^R , (66) 
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G 2 = 4( fet(e) Faf) £ a ^ S 9X-y 9p8 = V^9 R\p,p,u R Xp '^ ', 

Gs = T C a £ <T» e^ s e ABEF e CDGH e B ef e G x ef = y/=g R^ p , 

G, = \ tT* e ABCD T AB 7*jg = (R 2 - 4 R Xfl + R Xp ^ RT»» 
G 5 = T AB T° a f e^ s e A ef e c x e D p = y/=g (r 2 - 4 R Xfl R^ + R Xp ^ v B*#>>) , 

n _ J_ ( pua/3 ABCD -pAB C D\ ( uiua/3 -pAB A B\ _ \ppv p n 

7 ~ 4 e \ a & J V ^ a ^ / ~~ nrL \p,/J,is, 

/~i pva/3 -pAB -pAB fivafi p f?^ 3 

^8 — e pu cx/3 ^Xpj/ja/ 11 afii 

n — _ TT^B ttCD c /ifAp ^a/37<5 ABCD _ Xpj6 p p /if 

~~ 4det(e) pu al3 ® Xl ^ pS ~ n ^p,»v n -yS > 

^ _ 1 x-AB F CD pu\p afSjS CDEF A B E F _ Xpa/3 R r>pv 

~ det(e) ^ "J 8 7 5 A ~~ n \p,p,v n a/3 ■ 

Invariants 7-10 are P, T-odd, the rest are even. and G§ are full derivatives even if torsion 
is non-zero. The P, T-even _ S G» have been first c_ed by Nevifie 0. 

In the zero-torsion limit one replaces R K \^ V R K x,pu which satisfies the relations ( )44l) . 
Therefore, in this limit one has G 2 = G 3 , G 4 = G 5 , G 6 7 = 0, G 8 = G 9 = G w . Thus, in 
the zero-torsion limit one is left, apart from two full derivatives, with only two well-known 
invariants, namely 

V=g~R 2 and ^R kX R kX . (67) 



C. Invariants linear in torsion 

Using the covariant derivative of the curvature it is possible to construct invariants that 
are linear in torsion and linear in VP. The general structure from which all invariants of 
this kind can be derived is 

j^ab[cd] [ef] [gh] = det(e) ( v A p a/3/y5 ) e BX e°» e Dv e Ea e F ? e G ^ e m (68) 
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where V is the covariant derivative with the no-torsion Christoffel symbol (I37| . It belongs 
to the 6 3 <S> 4 2 representation of the Lorentz group, which contains 20 singlets: 

L[ = det(e)T A e A V x R Xp ^ } L> 2 = det(e) T A ef V x R^ Xp , 

L' 3 = det(e) T A V e A V^R^ P , L' 4 = det(e) T A V e A V^R pu \ 

L' 5 = det(e)T A e A ^V u] R, L' 6 = det(e) T A e A V x R x[p ' u]p , 

L' 7 = e x ^T A e A V x R, L' 8 = e^ s T A e A V X R X %, 

L' 9 = e^ s T A e A p V A i? 7 / p , L' 10 = T A e Aa e^ s V x R a ^ s , 

L'u = T A e Aa gfo e^ s V X R^ , L' 12 = det(e) T A e AX V X RT, 

L i3 = d^Je) eXpali T A U e A Va-Rq/3, 7 5, L' 1A = fct^t^t 01 ^ 5 T A e A V^p^s, 

L 'i 5 = l^f XpjS ^ T A e A V x R a ^ 5l L' M = ^ K e A ^ ^R a ^ 5 , 

L' 17 = T A e A p V x R a f, L' 18 = e*"* T A e A p V X R^ 

L' 19 = T A e Ax g a ^e^ s V x R a ^ s , L' 20 = T A e Aa V^R a ^ s . 

However, many of these invariants are zero or reduce to one another when one takes into 
account the additional symmetries of the standard Riemann tensor R, 

Ra/3,^u = Rfiu,a/3, (69) 

e p R a p,iJ,v = or R a p,fj,v + R^fi,au + RiMxfiv = 0, (70) 
as well as the Bianchi identity, 

£ VA^aft^ = or ^ \R a p,fiv + ^aRpX^u + V ' pR Xa ,nv = 0. (71) 



Contracting Eq. ( J7TI) with the metric tensor one obtains the identities for the Ricci tensor 
and the curvature: 

V^iT^ = -2V [ vR u]a , d p R = 2V U R^. (72) 

It is also helpful to keep in mind that the covariant derivative of the metric tensor and of 
the combination (1/ det(e)) e mv are zero. 

We immediately find that L[ = L' 2 , L' 3 = L 4 , L 8 = L' g , L' 10 = L' n , L' 17 = L' 18 and L' 8 = 

because of Eq. ( )69|) . L' 12 = L' 14 = L' 15 = because of Eq. (!70|) . and L' g = L' 10 = L' n = L' 13 = 

because of Eq. ( TfTj) . The invariants L 1; L' 2 and L' 6 are proportional owing to Eq. fl72|) . 

The invariants L' 10 = L' n = owing to Eqs. (17QI71jl . Therefore, actually only four linear 
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independent invariants are left: 

L'{ = det(e) T* e A ^ d^R = - det(e) v x d x R, 

L » = e \ P „v t a e A = 4det(e) a x d x R, 

LI = det(e) V^R^ P = § det(e) v x d x R+\ det(e) t P)Py V^', 

L'i = e A ^ T* e A ° V x R ptT = § det(e) a x d x R + § e A ^ V x R pa , 

and no linear combination of these invariants is a full derivative. They can be recombined 
in a simpler way: 

Li = det(e) v x d x R, L 2 = det(e) a x d x R, 
U = det(e) t p ^ V^R^, U = t x ^ t p J V x R pa . 

Since these invariants are linear in torsion, they are potential sources of torsion even in 
the absence of fermions, including the reduced torsion part, t A . 

We end up the derivative expansion here. We do not consider four derivative terms of 
the type (VT) 2 , (VT)T 2 , RT 2 and T 4 (there are many dozens of such terms) since they 
lead to even smaller corrections to the Einstein equation than the T 2 terms listed in Section 
VIA and considered in the next Section, see the estimate in Section II. However, all four 
derivative terms are, by dimension (which is 4), on equal footing from the point of view 
of the ultraviolet renormalization of the theory about curved background with generally 
nonzero torsion. Therefore, they should all be included, for example, in the "asymptotic 



safety" approach 



10 



32|. 



In the logic of the effective Lagrangians, which we assume in this paper, the gravitational 
action is an infinite series in d 2 /Mp such that it makes no sense in studying the stability 
against small runaway fluctuations from flat space-time in the concrete p 4 order since the 
inverse propagator of the fields is an infinite series in p 2 /Mp. What makes certain sense, is 
to study the stability of flat space-time at vanishing momenta but that is decided by the 
two-derivative terms. In addition to the usual condition that the Newton constant (or Mp) 
is positive, the new requirement is that the eigenvalues of the T 2 matrix are positive: this 
condition is summarized in the inequalities (163 p . 

It should be added that in the Einstein-Cartan formulation, none of the thinkable 
diffeomorphism- and local Lorentz-invariant action terms is strictly speaking stable under 



large nonperturbative fluctuations of the frame and spin connection fields 



. This obser- 



vation strengthens the argument that the present-day gravitation theory is but an effective 
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low-energy one, and therefore it makes sense to study systematically order by order the 
possible effects of the higher derivative terms. 



VII. INDUCED FOUR-FERMION INTERACTION 



If we ignore the TV-R terms ( 1731) that lead to high-order terms (V-R) 2 if we exclude the 
torsion, we are left with terms quadratic in torsion fl62|) and terms linear in torsion coupled 
to fermion currents (I55I) . It is important that in the leading order only the a M ,t> M part of 
the torsion couples to fermions. In the next order, however, when a derivative is added, the 
reduced torsion part may also couple to fermions [231 ] . 



In the leading order, if one integrates out the torsion the reduced torsion vanishes, 
whereas the Gaussian integral over a M ,t> M produces the 4-fermion interaction Lagrangian 



£* 4 



det(e) 



2 fM 2 \\- - M 4 

\ aavv av 



A B A B (gi 2 M 2 V - 2gtglMl + g+ 2 M 2 aa ) 

+ V B V B fa 2 Ml - 2g-g-Ml + g~ 2 M 2 a ) 

+ 2A B V B (gtgIM* v - (gfg^ + g x gt)M 2 av + gtg 2 M 2 aL 



(74) 



= V=g [h AA A B A + h vv V B V B + 2h AV A B V 

where A B = ^^ B ^^i is the axial and V B = ^7 B \1> is the vector current. The dimensionless 
constants gf 2 are defined in Eq. (|55j) and the masses M a>v are defined in Eq. (1621 . The A ■ V 
interaction term is C, P-odd and T-even. 

Certain particular cases of this Lagrangian have been considered before. For example, to 
compare it with the paper by Freidel et al. [7] we take gf = — 1, = a (see Section IV. C) 
and the 'minimal model' values of the torsion masses f|64"|) with the identification l ~ — 1/7, 
Ml = 2/G. In this case our general Eq. (J7J]) reduces to Eq. (23) of Ref. 3- 



VIII. STRESS-ENERGY TENSOR FROM FOUR-FERMION INTERACTION 

If the a^,v^ masses are of the order of the Planck mass as in Eq. (1641) the 4-fermion 
Lagrangian ( 1741) leads to a correction to the cosmological equation of the order of p 2 /Mp 
where p is the characteristic momentum of the fermion matter, for example temperature. 
Therefore, it is a tiny correction unless p approaches Mp but then one has to take into 
account higher terms in the derivative expansion, that are being neglected. 
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As discussed in Section V, the addition of torsion to the General Relativity in the fermion 
sector enlarges the gauge symmetry of gravity from the Lorentz SU{2)\ J x S77(2)r group 
to the ?7(2)l x ?7(2)r group which we know must be broken by the spontaneous or explicit 
masses of the a M , vector bosons. However, these masses need not be of the order of the 



Planck mass but could be much smaller, say, of the order of 10 TeV 26] . In this case 
the 4-fermion interaction ( IT4")) could be an important correction in the epoch preceding the 
electroweak phase transition. 

Anyway, there is an interesting problem of evaluating the contribution of the 4-fermion 
interaction to the stress-energy tensor in the r.h.s. of the Einstein-Friedman cosmological 
equation. This problem has been addressed e.g. in Refs. llMl4| using the ideas of a "spin 



fluid" 



34J-|36| . We think that this approach is unsatisfactory. Particles with spin one-half 



are always quantum, for example there are exchange effects, and that cannot be mimicked 
by any semi-classical model. At some point in the above references one has to average the 
spin-squared operator < s 2 >. This quantity is replaced by 1/4, why not 3/4? 

Meanwhile, averaging 4-fermion operators over a fermion medium is a common problem 
in Quantum Field Theory. With other fermion interactions (temporally) switched off, the 
contribution of the 4-fermion average to the Lagrangian is given by two terms - the 'direct' 
(Hartree) term and the 'exchange' (Fock) term, corresponding to two possible contractions 
of the ^ operators by the fermion propagator G(p), see Fig. 1: 



< (*r^)(w 2 *) > = \j 



Tr(G(pi)r! 



d 4 p 2 



Tr(G(p 2 )r 2 ) 



(2tt) 4 z 
Ti(G( Pl )T 1 G(p 2 )T 2 



(75) 



(2n)HJ (2vr) 4 i 

where Ti^ can be arbitrary Dirac and fermion 'flavor' matrices. If other fermion interactions 
need to be taken into account, one has to 'dress' the propagators and the 4- vertex. 





FIG. 1: Two contributions to the average of 4-fermion interaction: Hartree (left) and Fock (right). 



We emphasize that averaging over the medium should be performed in the Lagrangian; 
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the corresponding correction to the stress-energy tensor is then obtained by varying the 
Lagrangian with respect to g^ u . Averaging in the equation of motion makes no sense. 

In what follows we illustrate the use of Eq. ( 175]) by taking a medium of non-interacting 
fermions. 



A. Averaging 4-fermion interaction in a non-interacting medium 

We consider one species of fermions with mass m at temperature T and chemical potential 
p which corresponds to certain charge density p, to be specified below. It should be stressed 
that in a relativistic theory there is no strict way of separating particles from antiparticles: 
In a heat bath both particles from the upper continuum and holes (antiparticles) from the 
lower continuum are excited; the chemical potential regulates the difference between the 
number of particles and antiparticles, which is the only well-defined quantity. 

Neglecting interactions one writes the free fermion propagator 

G(p) = 1 = - = m+ (f i + i ^nYf + PiY , 7g s 

m — p 1 m — (p + iUn)^ — Pi7* m 2 + p 2 — (p + iu n ) 2 

where u n = 2nT(n + |) are the (imaginary) Matsubara frequencies. Integration over po 
becomes a summation over Matsubara frequencies: 

J toi'" .x" 
The charge density of the fermion gas is given by 

d 3 p 



= 4TV [ ^ * + lUn = 2 [^-(— —) (77) 

where e = y/p 2 + m 2 . This is nothing but the difference between Fermi-Dirac distributions 
for particles (positive p) and antiparticles (negative p). The integral can be easily evaluated 
in certain limiting cases; in particular, in the ultra-relativistic case m <C p,T one has 

The second term dominates at high densities when the massless fermion gas becomes de- 
generate. One can extract the chemical potential as function of p from this equation. 
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The thermodynamic potential of the fermion gas is [3 

d 3 p 



n = - 2l/T O og [( e ^ +i )(" +i 



(20 

'7* 2 -,4 . !-rf..J 1 



- V i8o T+ 5 TV+ i^' i r +o(m) - (79) 

We now compute the Hartree (direct) part of the 4-fermion interaction keeping in mind 
that in our case Ti.2 = (7 M ,7 M 7 5 ), see Eq. (1741) . We note immediately that the axial current 
does not contribute to the Hartree part as Tr(G(p)7 M 7 5 ) = 0. It could have contributed were 
the chemical potential different for left- and right-handed particles but we do not consider 
this possibility here. For the vector current we have two independent loop integrals (and 
traces) each of which is exactly of the same type as for the calculation of the charge density. 
Therefore, we obtain 

< V B V B > Hart = (^^Hart = \p\ (80) 

< A B A B > Hart = 0, 

< A B V B > Hart = 0. 

In the Fock (exchange) part the loop integrals and sums over Matsubara frequencies again 
factorize, however the trace does not. We obtain 

< V B V B > Fock = \p 2 - (81) 

< A B A > Fock = -p 2 + —a 2 , 

< A B V B > Fock = 0, 

where we have denoted 

d 3 p 1 r d 3 p 



a 



4T£ 



(2vr) 3 e 2 - (n + iu n ) 2 J e{2Tif 
At m C /i, T one gets 



1 - 



e t + l e t + l 



(82) 



T 2 /i 2 

a = 



6 2tt 2 ' 

The a 2 terms combine with the 0(m 2 ) corrections to the charge density (1781) . 
The full, Hartree plus Fock, contributions are 

< V B V B >Hart+Fock = 2 , (83) 

a a B 2 r fFt 2 

< A B A >Hart+Fock = + i 

< A B V B >Hart+Fock = 0. 
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These expressions should be put in the 4-fermion Lagrangian (1711) and we get: 



C 



4— ferm 



-g 



h 



w - A P 



m 



a 2 \+h 



AA 



4 



P 2 + 



m 



-a 



(84) 



We stress again that p here is the charge density, i.e. the density of particles minus the 
density of antiparticles, therefore it is expected to be small in the cosmo logical context. 



B. Derivation of the stress-energy tensor from the Lagrangian 

If the Lagrangian is known, the corresponding stress-energy tensor if found by the general 
rule 

1 dC 

Q U u = 2 ; — + derivative terms. (85) 

The problem therefore is to establish the dependence of ( 1841) on the metric tensor. This can 
be easily done if one realizes that T and \x are actually the zero components of 4-vectors, as 
seen from their use in Eq. (176 p . 

Let us illustrate how this logic works by finding the correct stress-energy tensor in the 
leading term in the fermion action. We take the first term in the thermodynamic potential 
(179)1 ; the simple case of the Stefan-Boltzmann law for massless fermions, Q = — 77r 2 T 4 /180. 
The partition function is 



Z = exp f-y) = ex P U J dt J d 3 x(-tt/V) 



where t is now the Minkowski time. It means that the corresponding Minkowski Lagrangian 
C = —Q/V generalized to the case of an arbitrary metric is 

7tt 2 

C = 180 ^ ( T » T "9n 2 , T, = (T, 0, 0, 0). 
Using the general Eq. ( 185)) we find immediately that in the co-moving frame 
7vr 2 d/ A e e N 7tt 2 \ T 4 , fi = v=0, 



60 I i T 4 ; il = v = 12)3 



3 

which gives, of course, the correct energy density e = O o and pressure p = Q n = 6 22 = 
G33 = I for the ultra-relativistic fermion gas. In fact e = 3p holds true for any relation 
between T and [i in that gas, since the dilatational current is conserved when there are no 



dimensional world constants, hence = 0. 
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We now find the stress-energy tensor following from the 4-fermion Lagrangian (184"]) . It 
will not satisfy the relation e = 3p anymore since the couplings hyy and h AA are not 
dimensionless but are of the order of 1/M 2 where M is the mass of a^, bosons. Neglecting 
the fermion mass we have 



£ 4 " ferm = ( - A hw + ih AA ) V=S P'_=^ + ^6r + £z- (86) 



3, 1 



4 " 4 











m=0 ' 


P 2 



p 2 T A 2p A T 2 p 6 



m =o 9 9vr 2 9vr 4 



At first glance, when we promote p, T to be 4- vectors there is an ambiguity: one can write 
p 2 T A as (p • p){T ■ T) 2 or as (p • T) 2 (T • T) or their combination. However, the result is 
independent of the decoding as long as p and T remain parallel, which is the case in the 
co-moving frame. What counts, is the total power of the polynomial in p, T, which is 6 in 
this case. Using the general rule ( )85|) we obtain for ultra-relativistic fermions: 



4— ferm 



3 1 

hvv + -rh A A ) p 2 n {-9tw + 6 5ya 5 vQ ) (87) 

1 ' m=0 



m=0 \ 4 4 

implying e 4 - ferm = 5p 4_ferm . This equation of state for the 4-fermion piece can be inde- 
pende ntly checked by using the general thermodynamic relations, see e.g. Eq. (1.4) of the 
book [381]. Indeed, if Z = exp(— f2(/x)/T) is the partition function one has 

pressure p = 1 — , charge Q = 1 — - — , charge density p = — , 

oV dfx V 

0(T In 2) E 

entropy S = — , energy E = —pV + TS + uQ, energy density e = — . 

oT V 



In our case we have from Eq. ( 186|) 

'3, 1 



lnZ 4 - fcrm = (-h vv + -h AA ^ p' 



V 

m=0 T' 



From the above general relations one immediately finds 

3, 1 



4 4 



e 4- t erm (/i; T) = 5 I ^ hyy + 1^ ) ^ 5p 4 - fem (/i, T) 



m=0 



confirming Eq. (1H7|) . The equation implies that [i and T are used as independent variables. 

However one may wish to express the stress-energy tensor in terms of the charge density 
p. To that end, one has to solve the equation p = (T/V)(d\nZ/dp) with respect to p and 
to substitute the function p(p) into the stress-energy tensor. In general, In Z is a sum of 
the main piece (179")) and the 4-fermion piece (1551) . such that p(p) is a complicated function. 
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But if the 4-fermion piece is a small perturbation, one can easily find the linear correction 
to the zero-order p(p) following from Eq. f!78|) . In this case one obtains 



P 4 - i ™( P ,T) = -(h vv + 1 -h ^ d{p2y) 



4— 1 A - LAA J dV 
The relation e 4 ~ fcrm = p4-ferm -^ as Deen used in Refs. 



a 4 

Q=const. 



7^ + 7^)p 2 = e 4 - ferm (p,T). 
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14j as following from the "spin 



fluid" approach. We see, however, that it is valid only if the a term (181]) from the Fock 
exchange contribution is neglected, and if the 4-fermion term is a small perturbation to the 
main part of the stress-energy tensor. If the 4-fermion interaction is the leading term, as 
assumed for the early cosmological evolution in the above references, the equation of state 
becomes e = 5p. 

C. Can the 4-fermion interaction be observable? 

On the whole, we come to rather pessimistic conclusions with regard to the observability 
of the possible 4-fermion interaction induced by integrating out torsion. If the 4-fermion 
constants hyv-, hAA are of the order of 1/Mp as assumed in the gravitational approach like 
in the 'minimal model' discussed (and criticized) above, it seems hopeless since it becomes 
significant only at particle momenta p ~ Mp but at these momenta we do not know the 
theory at all and anyway the derivative expansion fails. In addition, if for some reasons the 
4-fermion interaction becomes large, it must be included into the equation of state and not 
treated as a perturbation 12j-[l4| when it is overwhelming. 

We have mentioned that the masses M of a M ,t> M bosons could be of non-gravitational 
origin and therefore be, say, of the order of 10 TeV. That would increase the torsion- 
induced 4-fermion interaction by 30 orders of magnitude as compared to the previous case. 
Nevertheless, it is still hardly observable. In an ultra-relativistic medium the correction 
of the 4-fermion interaction to the stress-energy tensor is of the order of (/i 2 /M 2 )T 4 , as 
compared to the main contribution ~ T 4 , see Eq. (|86|) . A chemical potential in the TeV 
range is hardly imaginable. 

In general, it makes sense to introduce the chemical potential only for conserved quantum 
numbers. There are many conserved quantities in the late epoch, such as quark flavors and 
baryon (B) and lepton (L) numbers. However, in a late epoch the density is small and 
the torsion effects are probably negligible. The earlier we go into the evolution the fewer 
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quantum numbers are conserved. It was thought some time ago that electroweak interactions 
break B + L but preserve B — L quantum numbers, but today it is believed that both are 
broken since the Majorana type of neutrino is preferable. It looks like there are no conserved 



numbers at all in the epoch preceding the electroweak phase transition [39[, save the electric 
charge for which the chemical potential is zero. Therefore, it may well happen that in that 
early epoch the only possible contribution to the 4-fermion interaction is due not to charge 
density p (which is zero) but to the Fock 'exchange' part having the m 2 piece, where m 
is the fermion mass, see Eq. (I8ip . It may become large if there are super- heavy fermions 
but then their contribution is suppressed by the Boltzmann factor exp(— m/T), unless the 
temperature is of the same order of magnitude. 

Finally, we should mention the possibility that there is no thermal equilibrium in the 
epoch preceding the electroweak transition, meaning that temperature is not an adequate 
quantity. The Matsubara propagator ( I76|) is then irrelevant and should be replaced by 

G(p) ' 



where r is the relaxation time. In this case 1/r replaces, qualitatively, T and fi in the above 
equations for the estimate of the average 4-fermion interaction which, in principle, can then 
become sizable. We also mention an interesting possibility that an interplay of the evolution 
out of thermal equilibrium, and the potential C, P violation by torsion may lead to the 
baryon asymmetry of the Universe. 



IX. CONCLUSIONS 



We have systematically listed all possible invariants that may arise in General Relativity 
when one includes torsion, following the guiding principle of the derivative expansion. These 
include all possible invariants quadratic in torsion (5 invariants), quadratic in curvature (10), 
and linear in torsion and linear in the covariant derivatives of curvature (4). In the fermion 
sector, we have derived four possible invariants with torsion coupled to the bilinear fermion 
currents. We do not limit ourselves to P, T-even invariants. Some of the invariants are new, 
although most of them have been considered by different people before. 

In the leading one-derivative order, only 8 components of torsion (out of the general 24) 
couple to fermions, which can be cast into the Abelian axial (a^) and vector (v^) fields. If 
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parity is conserved, the interaction of a^v^ fields with fermions possesses gauge U(l)i, x 
C/(1)r symmetry, in addition to the Lorentz gauge symmetry SU(2) L x SU(2) R . Linear 
combinations of a^,v^ are the gauge bosons of this additional symmetry. 

However, the bosonic torsion-squared invariants break explicitly this symmetry as they 
provide masses to the Abelian bosons a M , or their linear combinations. From this point of 
view, such mass terms may look unnatural: spontaneous breaking could be more aesthetic. 
Certain justification for writing terms quadratic in torsion comes from the fact that they 
appear anyway from expanding the Einstein-Cartan action. This is called the 'minimal 
model' for torsion. However, we have shown that the minimal model leads to a non-positive 
mass matrix for the a^v^ bosons, therefore the minimal model cannot be complete. 

Assuming on purely phenomenological grounds that there is a positive mass matrix and 
neglecting higher-derivative invariants, we integrate out the torsion field and obtain the 
effective four-fermion action. It contains, generally speaking, axial-axial, axial-vector and 
vector- vector interactions. The effect of the first one has been studied in the past, with regard 
to its application to the Einstein-Friedman cosmological equation, using the so-called "spin 
fluid" approach. We find this approach unsound since particles with spin one-half are always 
quantum (for example there are exchange effects) and that cannot be mimicked by any semi- 
classical model. We present a systematic quantum field-theoretic method to average the 
4-fermion interaction over the fermion medium, and perform the explicit averaging in the 
case of free fermions with given chemical potential and temperature. The result is essentially 
different from that of the "spin fluid" approach. 

We arrive to rather pessimistic conclusions concerning the possibility to observe any 
effects of the torsion-induced 4-fermion interaction. However under certain circumstances it 
may have cosmological consequences, see Section VIII. C, but this has not been worked out. 
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